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relation with the evolution of accretion disks * 
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Abstract. We investigate the evolution of X-ray stellar activity from the age of the youngest known star forming 
regions (SFR), < IMyr, to about 100 Myr, i.e. the zero age main sequence (ZAMS) for a IMq star. We consider 
five SFR of varying age (p Ophiuchi, the Orion Nebula Cluster, NGC 2264, Chamaeleon I, and r] Chamaeleontis) 
and two young clusters (the Pleiades and NGC 2516). Optical and X-ray data for these regions are retrieved 
both from archival observations and recent literature, and reanalyzed here in a consistent manner so to minimize 
systematic differences in the results. 

We study trends of Lx and Lx/iboi as a function of stellar mass and association age. For low mass stars 
(Af < IMq) we observe an increase in I/x/iboi in the first 3-4 Myr and a subsequent leveling off at the saturation 
level (Lx/iboi ~ —3). Slowly evolving very low mass stars then retain saturated levels down to the oldest ages 
here considered, while for higher mass stars activity begins to decline at some age after ~ lO'^ years. 
We find our data consistent with the following tentative picture: low mass PMS stars with no circumstellar 
accretion disk have saturated activity, consistently with the activity-Rossby number relation derived for MS stars. 
Accretion and/or the presence of disks somehow lowers the observed activity levels; disk dissipation and/or the 
decrease of mass accretion rate in the first few Myrs of PMS evolution is therefore responsible for the observed 
increase of Lx/iboi with time. 
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1. Introduction 

Since the advent of sensitive X-ray imaging observations, 
late type (>F), low mass stars are known to be sources 
of X-rays with luminosities spanning more than 4 orders 
of magni tude, from Ly = 10' ^^ or less to more than 10'^^ 
ergs/sec (jVaiana et al.l Il98l|) . Given the similarity with 
the solar emission (e.g. thermal emission-line X-ray spec- 
trum, flares) it is believed that X-ray emission from low 
mass stars originates from stellar coronae that are es- 
sentially similar to the so lar one {Lx,q = 0.2-5.0 • 10^^ 
ergs/sec. |Peres et al.ll2000|) . albeit sometimes much more 
active. 

In order to observationally investigate the physical ori- 
gin of the large observed spread in Lx, homogeneous stel- 
lar samples have been successfully employed to restrict 
the number of free parameters. Open clusters and Star 
Forming Regions (SFR) are ideal targets in this respect, as 



open clusters and associations: individual: p Ophiuchi, 



all stars have approximatively the same age and metallic- 
ity. The comparison between different regions then ideally 
permits to investigate the age evolution and the metallic- 
ity dependence of X-ray activity. 

For late-type slowly rotating stars of nearby field pop- 
ulation, coronal and chromosp heric activity ind i cators 



correlate with stellar rotation jPallavic ini et alj. 
and with convective zone properties l|Nove'r'etal] 
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IPizzolato et al.ir2003l and references therein). Such corre- 
lations however do not hold for fast rotating stars, whose 
observed activity levels appear to have reached saturation 
(e.g. Lx/Lhoi ~ 10"'^). X-ray activity has also been found 
to decrease with increasing stellar age as stars loose an- 
gular momentum through an activity driven magnetized 
wind. Such observed relations are the main evidence that 
activity in the main sequence (MS) is produced by a stel- 
lar dynamo whose efficiency depends, at least for slow ro- 
tations, on stellar rotation and convection (the so called 
a — Lo dynamo) . 

While MS activity can be convincingly ascribed to the 
stellar dynamo, the mechanism responsible for coronal ac- 
tivity in low mass, pre-main-sequcnce (PMS) stars is not 
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understood. During the last 20 years sever al close-by SFRs 
have been observed in the X-rays (e.g. iFeigelson et alJ 
ll993HDamiani fc Micelalll995HFlaccomio et al.ll2003 a"bY 
finding bright X-ray emission from most PMS stars. No 
conclusive evidence for the action of a MS-type dynamo 
mechanism has been found for these stars. Although the 
presence of a convective layer does seem to be important 
for the existence of an active corona (massive fully radia- 
tive PMS stars do not show coronal emission) , no relation 
has been generally found between the X-ray emission level 
of PMS stars and their rotation periods^ and/or convec- 
tive region depth (or turn-over times). Instead, a relation 
has been ofte n found with st ellar mass or bolometric lu- 
minosi ty (e.g. IFeigelson et aLl fl993. 200^ iFlaccomio et al.l 
parameters unrelated to the a — uj dynamo. It is 
therefore not clear at present whether a new mechanism 
should be invoked to explain coronal activity in the PMS 
(e.g. a different dynamo), or the same dynamo acting on 
the MS works here in a saturation regime, and/or addi- 
tional phenomena occur in the PMS coronae that hide the 
underlying correlations with the a — lu dynamo parame- 
ters. 

The picture is probably further complicated by mat- 
ter accretion phenomena occurring on many PMS stars 
through a circumstellar disk: either accretion or the pres- 
ence of the disk may play a role in determining the ob- 
served X-ray activity levels. Several studies have indeed 
found that accreting stars have lower mean X-ray ac- 
tivity levels, and more time-vari able emission respect to 
their non-accreti n g cou nt erparts IINeiihauser et all 11995 ; 



F]a,ccomio et all \2()0^\ici. IStelzer fc Neuhauser 



20011) . However, this point is controversial: other au- 



thors have failed to detect a significant difference in 
activ ity levels of a ccreting and non accreting stars 



ig stars 

(e.g. IFeigelson et all ,2002: Prci bisch fc Zinneckeil l2002t 
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IFeigelso n et al.l Il993l) . As discussed by iFlaccomio et all 
((2003q1 . we believe that the discrepancy in the results is 
due, rather than to differences in the data themselves, to 
the different sensitivity of the techniques used to com- 
pare the activity levels of the two classes. In this work 
we will therefore assume that accretion disks do influ- 
ence activity. The physical nature of such influence is 
presently not unde rstood. Som e hypotheses, propose d by 
iDamiani fc Micelal (|I995) and IFlaccomio et alJ l|200 3bl. 
speculate that accretion (and/or the associated outflow) 
may modify the geometry of the magnetic field that con- 
fines the coronal plasma, either reducing the fraction of 
the stellar surface available for closed magnetic loops, or 
modifying more radically the magnetic field, maybe with 
the involvement of t he magnetized inner pa rts of the cir- 
cumstellar disk (e.g. lMontmerle et alJl200fi() . 



^ An exception is the Taurus region (e.g. 
IStelzer fc Neuhause J liool . However it is also possible that 
the correlation between activity and rotation is in this case the 
consequence of the possibly more fu ndamental correlations be- 
tween activity and mass-accretion ( Flaccomio ot al. , 2003b c) 
and between rotation and accretion(e.g. .Bouvier et al...l997.V 



While the study of individual SFRs and open clus- 
ters allows us to study the correlation (or lack of corre- 
lation) of X-ray activity with parameters such as stellar 
rotation, mass and circumstellar accretion indicators, im- 
portant clues to the mechanism responsible for activity are 
expected to come from the comparison of SFR at different 
evolutionary stages. 

In this work we assume, to first approximation, that 
the differences in activity between stars belonging to dif- 
ferent associations can be attributed exclusively to their 
different evolutionary stages. In other words, we assume 
that, in different associations, stars of a given mass are 
statistically indistinguishable at the time of their forma- 
tion (at least in regard to activity) and that their evo- 
lution is not influenced by different ambient conditions. 
These assumptions may be questioned because, for exam- 
ple: differences in initial angular momentum and stellar 
metallicity may affect X-ray activity; differences in envi- 
ronment, such as stellar density, can affect properties of 
the stellar population, such as the frequency of mult iple 
system s and the circumstellar disk lifetime ( Haisch ct al.L 
l2001bl) which may both influence observed X-ray activ- 
ity levels. Here, neglecting these possible complications, 
we will directly compare results obtained for stellar asso- 
ciations of varying age, in an attempt to investigate the 
effect of stellar evolution on activity. More precisely, we 
will select stars in definite mass ranges and we will follow 
the evolution with time of two activity indicators, Lx and 
Lx/Lhoi- 

In order to complete our program, we consider here 
six young PMS stellar groups with ages ranging from ^ 
1 to ~ 7 Myrs: p Ophiuchi, the Orion Nebula Cluster 
(ONC), split in the central and outer regions (see be- 
low), NGC 2264, the Chamaeleon I association and rj 
Chamaeleontis. In addition we consider two older ZAMS 
clusters: the Pleiades and NGC 2516. These regions have 
been selected so to probe a significant age range and be- 
cause optical and X-ray data suitable for our analysis were 
available either in the literature or in data archives. In all 
cases we try to consider stellar samples that are likely 
free from X-ray selection effects and, when needed, we re- 
analyze the published data in an homogeneous way. The 
latter is a crucial step because different instrumentation, 
assumptions and analysis techniques employed to derive 
physical quantities from both X-ray and optical observa- 
tions can result in significant systematic differences that 
would interfere with our purpose. 

In the following we first introduce the procedures em- 
ployed to consistently analyze (or re-analyze) the data 
for the various regions (Sect. I^J; we then discuss each of 
them individually, the sources of our data, the analysis 
performed, and the results obtained in regard to the de- 
pendence of activity on stellar mass (Sect. ISJ. Finally we 
present our evolutionary picture (Sect. ^ and propose a 
simple physical model that tries to explain our observa- 
tional findings (Sect.EJ. 
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2. Generals of data analysis 

The quantities we will deal with are two X-ray activity in- 
dicators {Lx and Lx/iboi), stellar masses and ages. The 
estimates of these quantities are often subject to signifi- 
cant uncertainties, due both to errors in the observational 
quantities and to the assumptions made in order to con- 
vert these latter into physical parameters. In order to make 
the comparisons of activity among different associations 
meaningful, special care must be therefore taken on the 
one hand to reduce random uncertainties and on the other 
to minimize systematic differences due to the non-uniform 
approaches and assumptions made in the different studies. 
In the following two subsections we review the steps taken 
to this effect. 

2.1. X-ray luminosities 

We utilize data obtained with four X-ray detectors: 
ROSATRRl and PSPC, Chandra HRC-I and ACIS-I. In 
the original works from which these data are taken. X-ray 
luminosities were calculated from measured instrumental 
count-rates assuming different spectral bands, the choice 
depending largely on the spectral response of the instru- 
ment. These X-ray luminosities are not directly compara- 
ble. The count-rate to luminosity conversion factors were 
computed in all cases assuming optically thin isothermal 
plasma emission. The assumed plasma temperatures (kT), 
however, varied from case to case, as well as the way in 
which interstellar extinction was accounted for. Here we 
recompute all conversion factors, using PIMMS^, for a 
single spectral band (0.1-4.0 keV). For PMS stars we as- 
sume kT — 2.16 keV, consistent with recent results on 
severa l SFRs (e.g. iFeieelson et al.l|2002l: iFlaccomio et all 
l2008al llmanishi et al.ll2001af) . and an hydrogen column 
density, TVh, proportional to optical/IR extinction, Ay or 
Aj, where the latter is measured, with the exception of rj 
Chamaeleontis (Sect. IX^ . individually for each member 
of the association. 

For the two MS clusters, NGC 2 516 and the Pleiades, 
we assume (cf. iGaene et al.l ll995b^ lower temperatures, 
between ^ 0.4 and ~ l.OkeV. Given the small and quite 
uniform interstellar extinction, A^h was assumed, for stars 
lacking individual measurements, to correspond to the 
mean Ay of the entire cluster. 

Departure from our assumptions for kT and A^h will 
result in errors in the count rate to Lx conversion, whose 
magnitudes depend on the detector used for the observa- 
tion and on the range of variation of the two parameters 
for the given stellar population. In the following we will 
discuss, for some of the regions considered, the sensitivity 
of the conversion factor to these assumptions. In all cases 
the likely systematic uncertainty in Lx is smaller than 
~ 0.1 dex. 



2.2. Masses and Ages 

We make use of ISiess et all tOOd SDF hereafter) evo- 
lutionary tracks'^ in order to estimates masses and ages 
of PMS stars from effective temperatures and bolomet- 
ric luminosities. These latter values are taken from pub- 
lished works, derived in most cases from medium resolu- 
tion optical or IR spectra. For the two oldest PMS regions, 
Chamaeleon I and 77 Chamaeleontis, estimates are based 
on optical/IR phot ometry. Ages for the two ZAMS clus- 
ters are taken from iMevnet et al.l lIlQQ,*^ and are derived 
from the turnover point in the upper main sequence; we 
then estimate stellar masses from optical photometry and 
the theoretical isochrones appropriate to the clusters age 
(SDF; solar metallicity) , transformed from the theoretical 
^boi - Toff plane into the ob servational V - (B — V) plan e 
with the relations given by iKenvon fc HartmannI l) 19951) . 
Note that masses for both PMS and MS stars are esti- 
mated from the same theoretical evolutionary model. In 
order to alleviate the effect of residual uncertainties in 
the mass estimates (e.g. due to the different temperature 
scales used in the original studies the data are taken from) 
we will consider wide mass range s , simila rly to what done 
for the ONC bv iFlaccomio et alJ ll2003bl) . 

We now turn to the problem of PMS stellar ages: in- 
dividual ages inferred from PMS evolutionary models are 
known to be subject to large uncertainties. In addition to 
the uncertainties in the theoretical evolutionary models, 
bolometric luminosities are subject to an art ificial spread 
due to a number of effects jHartmannl l200l(l . As a result 
the real dispersion of stellar ages in a given star form- 
ing region may be much smaller than estimated from its 
HR diagram. We decided to characterize each PMS stellar 
group with a single age, i.e. the median of the ages derived 
from the SDF tracks. As long as the shape of the age distri- 
bution due to the artificial spread in Lboi is at least similar 
for all star forming regions, we can assume that the median 
of the measured ages is in some close and unique relation 
with the real median age. Figure ^ shows the distribution 
of individual stellar ages, as derived in the following sec- 
tion, for each of the six PMS stellar groups considered 
here. Also indicated are the median values which will be 
used in the following. 



3. The star forming regions and open clusters 

We now discuss the eight stellar groups taken in consider- 
ation, describing the stellar samples and the analyses per- 
formed, along the lines sketched in the previous section, 
in order to make the results comparable to each other. We 
proceed in chronological order, starting from the youngest 
stellar association, according to the sequence suggested by 
Fig. [2 



^ Portable, Interactive, Multi-Mission Simulator 



^ See also the web site:,http://www-laog. obs.ujf-grenoble.fr/activites/starev 
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- p Ophiuchi 5.99 (43) 

- ONC Trapezium 5.17 (130) 
■ ONC outer 6.28 (545; 

NGC 2264 6.34 (l 78) 

- Chamaeleon I 6.60 (71 

- Tj Chomoeleontis 6.86 (10) 




Loq(Age) [y,.] 

Fig. 1. Cumulative age distribution for the six PMS re- 
gions under study. The legend on the upper-left side re- 
ports, for each region, the median of the logarithm of the 
age (in years) and, in parenthesis the number of stars in 
the sample. 



3.1. p Ophiuchi 

p Ophiuchi, with an age of < 10^ years, is the youngest re- 
gion here considered. It is a medium size association form- 
ing out of a dense molecular cloud, and most of its mem- 
bers are deeply embedded in the cloud. It has been there- 
fore studie d with more success in the near IR and X-ra y 
bands fe.g. lLuhman fc Rie ke"l999':'lmanishi et all2001a|) . 
Our source o f optical/IR data f or p O phiuchi members is 
the work by iLuhman fc who obtained NIR 

spectroscopy and photometry of 114 stars and derived, 
for 69 of these, effective temperatures, J band extinctions 
and bolometric luminosities. We restrict our analysis to 
the field of view (FOV) of the two Chandra X-ray ob- 
servations described below, so lowering the census of our 
well characterized members (i.e. with T^s and iboi) to 51. 
We estimated masses and ages for 43 of these stars using 
SDF evolutionary tracks, the remaining eight having ap- 
parently masses lower than the lowest mass for which SDF 
tracks are available (O.IM©). These 43 stars comprise our 
reference sample for p Ophiuchi. 



3.1.1. X-ray data and analysis 

We have retrieved, from the Chandra Data Archive 
f htt p : // asc . harvard . edu /cda/ ) , two deep (w 100 ksec 
each) and slightly overlapping ACIS-I observation of 
the p Ophiuchi cloud (Obs. Id. 635, 637). A thorough 
analysis of the first obse r vati on has been already re- 
ported by llmanishi et alJ l(2001ai) . focusing mainly on 



mediu m resolution X- r ay spe ctra and on source vari- 
ability; llmanishi et alJ l|2001bj) have used both observa- 
tions to study X-ray emission of the known brown dwarfs 
and candidate s brown dwarfs in the region. However, 
llmanishi et alJ (1200 Ihl) do not produce a list of sources 
and llmanishi et al.l l|200iaj) do not compute upper limits 
to X-ray luminosities of undetected cluster members, a 
necessary step in order to avoid selection effects. To uni- 
form our analysis to the general principles stated in Sect, 
n we re-analyzed the two dataset. 

First, in order to selectively lower the background level 
we filtered the two event files retaining only events with 
energies in the 0.3 to 8.0 keV range, likely containing the 
largest fraction of the X-ray flux from p Ophiuchi mem- 
bers. We then dete cted sources and meas ured count rates^ 
using PWDETECT ijDamiani et al.1 ll 9971 Damiani et al., 
in preparation)^ and a signal to noise ratio threshold of 
5.10, which, given the background level of the two ob- 
servations, is expected to yield ~ 1 spurious detection in 
each FOV. Ninety sources were detected in obs. 635 and 
77 in obs. 637, after removal from the latter source list 
of ~20 spurious detection due to a bright CCD column in 
coincidence with a lum inous X-ray source. For comparison 
llmanishi et alJ l)2001a(l detect 87 sources in obs. 635. Seven 
sources located in the small overlap region of the two ob- 
servations (~ 40arcmin^), were detected in both observa- 
tions. In the same region 6 sources appeared instead in 
only one of the two datasets. We created a merged source 
list, assigning to the sources detected in both observations 
uncertainty-weighted averages of count rates and spatial 
coordinates. Our merged source list comprises 160 distinct 
X-ray sources with count rates ranging from 6.3 • 10^'"' to 
0.79 cuts • s-^ 

We cross-identified our detected X-ray sources with the 
list of 43 selected p Ophiuchi members discussed in the 
previous section. For comparison purposes (see below), we 
have also identified, with the same list of membe rs, the X- 
ray sources detected bv llmanishi et al.l l|2001aj) . In both 
cases the identification process went smoothly, thanks to 
the high spatial resolution of Chandra and the relatively 
low spatial density of objects in the region. We identify 36 
members with 35 X-ray sources, a close pair of members 
(GY240A and GY240B) being not resolved in the X-ray 
image (we treat the measured count-rate of the pair as an 
upper limit to the count-rate of each star). For the seven 
members undetected in the ACIS-I data we computed up- 
per limits to count rates using pwdetect and the same 
signal to noise ratio threshold used for detections. 

We estimate X-ray luminosities of stars in our sample 
from ACIS-I count rates and assuming a distance of 165pc 



* Effective exposure times were taken from an expo- 
sure map computed using the tools in the Chandra 
Interactive Analysis of Observations (ciAO) package 
(http://asc.harvard.edu/ciao/), for an incident energy of 
2.16 keV. The choice of energy is however not critical as 
vignetting depends only mildly on energy for the range of 
energies relevant to stellar coronae. 

^ See also ,http: / / www.astropa.unipa.it / progetti_ricerca/PWDetect| 
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ijDame et al.lll987j) .^ Count rates to flux conversion fac- 
tors were computed, as described in Sect. 12.11 i.e. taking 
kT = 2.16keV and A^h proportional to the J band absorp- 
tion, Aj (see below). Thanks to the anal ysis of medium 
resolut ion ACIS-I spectra performed by llmanishi et alJ 
which also yielded independent estimates of X- 
ray luminosities, we are in the position to check these as- 
sumptions. 

The median(m ean) of the quiescent s tate kT derived 
from spectral fits ijimanishi et all l2001a|) for 28 stars in 
our member list is 2.1(2.3) keV, indeed similar to our 
value. Moreover, considering the relevant uncertainties of 
these best fit kT, their distribution is consistent with a 
narrow range of values: '^50% fall inside the kT — 1.5 — 2.5 
keV range and another 35% have 90% error bars consistent 
with the same interval. 

In regard to extinction. Fig. [21 shows the proportion- 
ality between X-ray derived A'h and Aj. The proportion- 
ality constant was assumed to be the median of the mea- 
sured A^hMj: 4.72 ■ 10^^. This is somewhat lower than 
what would be expecte d on the basi s of a normal opti- 
cal/IR extinction curve l|Mathid . lT990|) . Ay/Aj=3.55j^nd 
the re lation between A^h and Ay given either by iRvter I 
lll996t) . TVh = 2.0 • 10^1 Ay or IPredehl fc SchmittI mm . 
Nn = 1.79 ■ lO'^^Ay. With these assumptions we would 
have derived A^hMj (ArHMv)(AvMj) = 7.1 • lO^i and 
6.35 ■ 10^^ respectively, i.e. about 30-50% higher than our 
value. This would have in turn resulted in overestimating 
the conversion factors (and Lx) by 0.1-0.2 dex, depending 
on absorption in the Ay = 2-50 range. The departure from 
the normal relation may points toward a pe culiarity of 
the extinction law in the p Ophiuchi c loud fcf. IVrba et al.l 
ll993HKenvon. Lada. fc Barsonvlll998|l . 

Finally, Fig. El compares our X-ray lu minosities with 
those estimated bv llmanishi et all (2Q^^ on the basis of 
spectral fittings; 90% confidence interval are also shown 
for the latter. More th an one estimate of Lx is given by 
llmanishi et al. I ^l2001a^ for some variable sources corre- 
sponding to quiescent and flaring states (indicated with 
different symbols in the Figure) while we only derive an 
averaged Lx over the whole observations. In general , the 
estimates for quiescent Lx bv llmanishi et all l)2001a|) and 
our estimates are fully compatible, with the possible ex- 
ception of some flaring sources, for which the difference 
can be easily attributed to the presence of the flare. 

We conclude that, at least for p Ophiuchi stars, the 
assumption of a single plasma temperature for all sources 
and of optical extinctions to calculate conversion fac- 
tors yields realistic values of Lx- From the dispersion of 



® There is some uncertainty on this distance: an Hipparcos 
based estimate of the distance to the Upper Scorpius re- 
ion, of which the p Ophiuchi cloud is part, indicates ~145pc 
Ide Zeeuw et a l.. 1999'). Had we adopted this value Lx would 
be lowered by ~ 0.12dex. The Z/x/iboi values, on which our 
main conclusions are based would remain unaffected. Note also 
that the median stellar age of the region members would in- 
crease by ~0.17dex, becoming almost equal to that of the ONC 
Trapezium region (Sect. lT^ . 




-J 10 15 

Aj [r-ag] - Luhman & Rieke (1999) 

Fig. 2. A^H derived from fittings of X-ray spectra 
l lmanishi et all boOlal) vs. J band extinction. Aj, mea- 
sured for the same stars bv lLuhman fc Riek^ l|l999l) . The 

straight line corresponds to the median of the N-^/Aj ra- 
tios: 4.72 • 10^^ Error bars are 90% confidence intervals 
derived from spectral fits. 



spectral-fit kT and A'h about our assumed values, we es- 
timate that the mean error on our conversion factors is of 
the order of 0.1 dex. 



3.1.2. Activity vs. mass 

We study the relation between activity and stellar mass 
using the same method described in iFlaccomio et all 
|2003ll) for the Orion Nebula Cluster. Figure^ shows X- 
ray luminosity functions for the stars in five mass bins. 
The same Figure shows the median Lx and the dispersion 
(25 and 75% quantiles) of the luminosity functions as a 
function of mass. Figure El shows the same kind of anal- 
ysis for our second activity indicator: Lx/Lhoi- We draw 
the following conclusions: 1) the median Lx is correlated 
with mass; 2) the median Lx/iboi (10~'*'^-10~^ °) is sig- 
nificantly below the canonical saturation level, 10~^; 3) 
ix/iboi is correlated with mass, at least at for M < IMq. 
Note that the correla tion of Lx I i^hni with ma ss was also 
found in the ONC bv IFlaccomio et alJ l)2003b^ : here how- 
ever the dependence appears to be somewhat steeper. 



3.2. Orion Nebula Cluster 

The X-ray (Chan dra HRC) data for the ONC are de- 
scribed in detail in Flac comio et !n]i2003a). There we also 
defined an optically selected, extinction limited sample of 
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29.0 29,5 30-0 30-5 

Log(Lx) [ergs s"'] - this work 



Fig. 3. Comparison of X-ray luminosities estimated from 
the sam e ACIS-I data in this work and bv llmanishi et all 
ll2001a^. Filled circles indicate that the Lx estimate by 
llmanishi et~aLl l)2001a|) refer to the quiescent source emis- 
sion, while empty triangles refer to flares. The numbers 
be side each point i s the so urce identification number given 
bv llmanishi et all l)2001a|) . Note that up to three estimates 
(one quiescent and two flares) are given for each source. 
Error bars are 90% confidence intervals derived from spec- 
tral fits. 



well characterized ONC members (there referred to as the 
optical sample) that we also use here. 

For the present purpose, i.e. the study of the evolution 
of activity in the PMS, we split the region surveyed by 
the HRC observation in two sub- areas: the central 5' x 5', 
roughly identifiable with the so called Trapezium region, 
and the re maining of our HR C FOV. It has already been 
noted fe.g. lHillenbrandlll997l) that the Trapezium region 
is of somewhat more recent formation respect to the outer 
part of the molecular cloud. FigureHshows the cumulative 
age distribution for the two regions: it is evident that, 
although the two distributions overlap there is an excess 
of very young stars in the Trapezium. Part of the overlap 
may be due, other than to the artificial and real age spread 
in the two regions, to the fact that the area we identify 
with the Trapezium, likely contains, for a projection effect, 
a number of stars that are in fact located in front of the 
Trapezium and should belong more properly to the outer 
region. In any case, following the discussion in Sect.Elwe 
will assume that the median activity levels of stars in these 
two regions are representative of the two median stellar 
ages. 

Figure 21 and [5l show our standard analysis of activity 
vs. mass for the two spatially segregated regions. We note 
that the same qualitative trends found for the whole ONC 



by iFlaccomio et alJ l)2003bl) are retrieved in each of the 
two sub-regions. 

3.3. NGC 2264 

A recent analysis of stellar activity in NGC 2264 was pre- 
sented bv IFlaccomio et all ||2003 ct: thes e authors utilized 
published X- r ay and optical data ( Flaccomio et all l200(i 
iRebull et al.L l2002bl) to investigate the relation between 
circumstellar environment and activity. The data were 
there reanalyzed following the general criteria stated in 
Sect. 121 Here we use the same data and the same refer- 
ence stellar sample. The sample comprises 178 candidate 
membe rs from the spectroscopic sample o f £ebull et alJ 
l|2002bl) . to which iFlaccomio et alJ l(2008c() were able to 
assign masses, ages and X-ray luminosities (or upper lim- 
its), thes e latter derived from the ROSATYHM count rates 
listed bv IFlaccomio et all l)2000|) . 

Figure 21 and [51 show our analysis of Lx and -^^x/iboi 
vs. stellar mass. Given the large number of upper lim- 
its, it is quite clear that the sensitivity of the presently 
availa ble X-ray data {ROSAT from 'Flacco mio et alJ 
is not sufficient to define the low luminosity tails of 
the X-ray luminosity functions so that the median Lx is 
defined only in two mass bins and the trend of Lx with 
mass seen for the other regions is not so apparent here. 
The median Lx/Lboi shows a pattern similar to the other 
regions, being close to the saturation level for M < I.OMq 
and declining for larger masses. 



3.4. Chamaeleon I 

The Chamaeleon I association, with an age of 5-6 Myrs, 
is the second oldest star forming region considered in this 
study. Like for NGC 2264, the optical and X-ray data 
we use here are described in detail by IFlaccomio et alJ 
(|2003cj) . where the activity/disk connection was analyzed. 
The original data t hat work is based on were taken for 
the most part from iLawson et al . (1996), but reanalyzed 
consistently with what done for the other SFRs. 

Figure 21 shows the X-ray luminosity functions in our 
chosen mass ranges and the Lx vs. mass scatter plot. 
Figure refers, in the same format, to Lx/Lboi- The gen- 
eral trend of increasing Lx with increasing mass, seen in 
the other star forming regions, can be clearly observed 
here. Lx/Lboi seems to be close to the saturation level 
at all masses. In particular we note that, respect to the 
younger association, p Ophiuchi and the ONC, we do not 
observe a decrease of the median Lx/Lboi at very low 
masses. 



3.5. 7] Chamaeleontis 

rj Chamaeleontis is a small, nearby (97 pc) and relatively 
old (- lOMyr) PMS cl uster discovered t hroug h ROSAT 
X-ray observations by iMamai ek et a,lJ lll999l i. The X- 
ray properties of 12 members have been investigated by 
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iMamaiek et alJ l)2000() . while iLawson et all l)200l|) report 
optical photometry and derive rotational periods. The to- 
tal r] Chamaeleontis known population counts 15 stars: 3 
are early type stars (one of which, HD75505 is undetected 
in the X-rays), 11 are WTTS and 1 is a CTTS. Ten of 
the 12 low mass membe rs have been discovere d through 
detection in the X-rays l|Mamaiek et alJ.ll999|) . The fol- 
lo wing optical search fo r additional members conducted 
bv iLawson et al.1 l|2002|) resulted in the discovery of only 
two non X-ray selected members, among which the CTTS 
(ECHA J0843.3-7905). 

In the following, reassured by the high detection ef- 
ficiency among the known low-mass cluster population, 
we will only consider X-ray detected mem bers. The X-ray 

i ROS ATiniT) count rates are taken from lMamaiek et al.l 
200(1) while the optical data, Tes and Lboi, come from 
ILawson et al.l (|2001). Following the previous investiga- 
tions we assumed zero absorption along the line of sight. 
The determination of masses, ages and X-ray luminosities 
followed the principles stated in Sect. El 

Figure 01 shows the X-ray luminosity functions in our 
mass ranges and the trend of median Lx with mass. Figure 
prefers, in the same format, to Lx/iboi- We note that, 
notwithstanding the low number of stars, the data points 
are compatible with the same qualitative results found 
for low mass stars in the slightly younger Chamaeleon I 
region: a dependence of Lx on mass and a constant 'sat- 
urated' Lx/ibol- 

3.6. The Pleiades 



The Pleiades open cluster, with an age of ~ lOOMyr 
i Mevnet et a 

filial and a distance of only w 127pc, is 
possibly the best studied young ZAMS cluster in the so- 
lar neighborhood. Its stellar population has been nearly 



completely identified and several X-ray studies perfo rmed 
wiihROSAT fe.g.lStanffer etTillTflflilMicela et 3,1.1119961 

I1999I) have measured X-ray luminosities for a representa- 
tive fraction of Pleiades members. Here we use this wealth 
of information for the purpose of extending our study of 
the evolution of activity from the pre-main-sequence to 
the ZAMS. 

Our sample of Pleiades members includes ^ 500 stars 
listed in the open cluster database of C.F. Prosser and J.R. 
Stauffer'' and that fall in the FOV of the X-ray observa- 
tions described below. We estimate stellar masses from the 
theoretical SDF isochrone and Vq magnitudes as explained 
in Sect. 12.21 The choice of using the V magnitude instead 
of, for example, the B — V color as a mass indicator is due 
to the fact that V is available for all the stars in our sample 
(vs. ^ 54% for B — V). We exclude from the following anal- 
ysis 38 stars that, in the Vq vs. {B — V)q diagram, are more 
distant than 1 magnitude from the isochrone (along the 
Vb axis). We note that, given the non-perfect adherence of 
the observational member locus to the SDF isochrone and 
to the presence of unresolved multiple systems (having 



,http://cfa-www.harvard.edu/^stauffer/opencl/index.html| 



lower Vb than those pertaining to the primary alone) our 
mass estimates are somewhat uncertain. We get a rough 
idea of the uncertainties on mass estimates by compar- 
ing our masses derived from Vb with those obtained from 
the same isochrone through {B — V")o. The difference in 
the individual masses are typically smaller than 0.05 dex, 
small compared to the width of the mass ranges used to 
compute the Maximum Likelihood distributions of Lx and 

ix/ibol- 

Bolometric luminosities, used for computing Lx/iboi) 
are estimated, in the same way as masses, from the theo- 
retical isochrone. 

X-ray data are taken from the published analysis of 
16 different X-ray exposures performed with ROSAT, 8 
using the P SPC Micela et alJ }l998) and 8 using the HRI 
l)Micela et al.l Il99^ . Each star in our reference catalog 
was observed between one and nine times (~ 3 times on 
average). However, possibly due to sensitivity limits, only 
46% of the stars in our reference list are detected, while 
for the remaining ones only upper limits to the X-ray lu- 
minosities are available. 

The es timates of X-ray luminosities we use are 
taken from 'Micel a et aP l)l998|) for the PSPC data and 
iMicela et al.. (,19991) for ;he HRI data. Both refer to the 
0.1-2.5keV spectral band and are based on the following 
assumptions: a coronal temperature kT = 0.8keV and ab- 
sorption A^H = 2 ■ 10^1 • 3.1 • E{B - V), where E{B - V) is 
measured individually for about 15% of the stars, while for 
the remaining ones a mean value, 0.04, is used. The spec- 
tral band luminosities refer to is narrower than our stan- 
dard 0.1-4.0keV band; anyway, because of the low coro- 
nal temperature of Pleiades sources, the difference in con- 
version factor due to the difference in band is small: for 
kT — O.SkeV, the fraction of emitted energy falling in 
the 2.5-4.0keV band, and therefore the amount by which 
we underestimate Lx, is < 0.01 dex. Only for kT ^ 1.5 
keV, hotter than typ i cal MS coronal sources l|Gagne et all 
Il995hl: jPaniel et all l2003() . the difference in Lx due to 
the narrower spectral band reaches ~ 0.05 dex. Softer 
spectra, kT ~ 0.4 — 0.5 keV, have been inferred for main 
sequence F type stars llGagne et all Il995bt iFlaccomiol 
l2002t iDaniel et all l2003|) : assuming kT = O.SkeV results 
in this case in an overestimation of Lx of up to ~ 0.05dex, 
which is small compared to the other sources of random 
uncertainties (e.g. counting statistic and distance spread). 
The assumption of an average value of absorption, for stars 
for which no individual measurement of E{B^V) is avail- 
able, is also justified: 80% of the stars for which E{B — V) 
has been measured have a corresponding A^h ^ lO'^^ cm^^, 
which, for kT ~ 0.8 keV translates into a difference on the 
individual conversion factor of < O.ldex respect to assum- 
ing the mean absorption. Moreover, while stars that have 
individual extinction estimates are located in a sky area 
of high nebulosity (this being the reason that prompted 
the measurement of individual extinctions), the remaining 
stars are located in an area apparently devoid of clouds 
and most likely have fairly uniform E{B — V); the former 
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estiniate of 0.1 dex for the uncertainty on the conversion 
factor is therefore most hkely too pessimistic. 

Given the multipHcity of observations, in order to 
assign a single Lx value to each of our stars, we pro- 
ceeded as in Flacco mio et al.l l)2000|) . i.e. we compute max- 
imum likelihood mean luminosities using all the avail- 
able information (detections and upper limits). This was 
implemente d by computin g, with the ASURV package 
ijFteigelson fc NelsonlllQSfil) . the Kaplan Mayer estimator 
of the values of log(Lx) available for each source. Because 
of the time variable nature of X-ray emission, the thus 
derived mean log(Lx) of ~ 6% of the stars which were 
detected in some of the observations turned out to be up- 
per limits. Out of the 483 stars in our reference sample we 
therefore have estimates of mean log(Lx) for 40% of the 
stars and upper limits for the remaining ones. 

Figure shows X-ray luminosity functions in our cho- 
sen mass bins and the log(Lx) vs. mass scatter plot. The 
same analysis is repeated for log(Lx/iboi) in Fig. [7| Low 
mass stars (log(M/MQ) < —0.3) appear to be saturated, 
with log(ix/iboi) ^ ^3, and show the consequent mean 
correlation of Lx with mass. Higher mass stars however 
are not saturated, so that the relation of luminosity with 
mass is lost, as also evinced by the onset of the activity - 
rotation connection which is indeed observed at this mass 
boundary fcf. .Micela et al. . 1999). 

3.7. NGC2516 

The open cluster NGC 2516 is in many respects simi- 
lar to the Pleiades, and it has indeed been often named 
the so uthern Pleiade s . It is however slightly older, ~ 140 
Mvrs jMevnet et all Il993|) . and more distant (387 pc, 
I.Teffries et al.lll997l) . 

Our data are b ased on the photometric survey of 

I.Teffr ies et ahl and the an alysis of eight C handr a 

X-ray observations reported by iDamiani et al.1 l|2003|) . 
Our initial member list is the same as that used by 
Damiani et a l. ( 2003) and comprises 725 stars indicated 
as members by I.Teffries et al.l l)200lh and falling in the 
FOV of at least one of the Chandra observations. Stellar 
masses, and bolonietric luminosities, were derived analo- 
gously with what described for the Pleiades cluster from 
Vq magnitudes and a theoretical (SDF) isochrone, the only 
difference being in the choice of the isochrone appropriate 
to the larger age of NGC 2516. We excluded from our 
reference sample 52 stars that, in the Vq vs. {B — V)q dia- 
gram, lie farther away than 1 mag (along the Vq axis) from 
the 140 Myr isochrone and four stars that were found to 
be outside the My range covered by the isochrone. 

We adopt X-ray luminosities from [Pamiani et alJ 
l|200.3l) . These are averages of the luminosities estimated 
from detected source counts in all the analyzed observa- 
tions. Conversion factors between counts and luminosity 
in the 0.1-4.0keV band are there computed assuming mii- 
form extinction, A^h = 8.2 • 10^°, and, for bright sources, a 
two temperature spectral model (one of which fixed at 



Table 1. Optical and X-ray data for the stellar sample 

N. Region Mass log(Age) log(Lx) log(Lx/Lboi) 

[Mq] [yr.] [ergs s'^] 

(1) (2) (3) (4) (5) (6) 



Notes - Column (1): sequential reference number; column (2): 
identifier of the region the star belongs to. The table, 
comprising 2129 rows, is available in electronic form at: 
://www.astropa. unipa.it/~ettoref/Act-EvoLTl. datl 



0.26keV), constrained by observed hardness ratios. For 
weaker stellar sources and upper limits, the conversion fac- 
tor is taken as the median of the values derived for bright 
sources with the same spectral type. Damiani ct al. (2 0031) 
estimate that, given the range of variation of the conver- 
sion factors derived from hardness-ratios, the uncertainty 
in individual log(ix), introduced by an inaccurate spec- 
tral model, is at most ^ 0.15 dex. 

Figure jSj and show our standard analysis of activity, 
Lx and Lx/Lboi, as a function of mass. The same quali- 
tative trends observed for the Pleiades are retrieved here. 

4. Evolution of coronal activity 

We now utilize the data presented in the previous section 
to study the evolution of stellar activity in specific stellar 
mass ranges. Tabled available in electronic format, sum- 
marizes the data described in Sect. O and also utilized in 
the following. For our purpose (see Sect. I2.2|l we assume 
that the median activity level of each roughly coeval PMS 
stellar group is representative of the PMS activity at the 
median of stellar ages in the same group. 

Figure |S| shows, in the left hand column, for each of 
our chosen mass ranges, the median log(Lx) of each of the 
eight regions studied in the previous section as a function 
of median association age. The same analysis is presented 
for log(-Lx/-Lboi) in the right hand column. In order to 
assess the significance of the difference between Lx and 
Lx/ Lhoi distributions of stars in a given mass range and 
belonging to different regions (i.e. having diff'erent ages) 
we have performed two population statistical tests using 
the ASURV package. ASURV performs five different tests 
that give quantitatively diff'erent results. In comparing dis- 
tribution of activity indicators we will consider, quite ar- 
bitrarily, significantly and possibly different two samples 
for which the majority of the tests give lower than, re- 
spectively, 1% and 10% probability that the samples are 
not drawn from the same parent population. In order to 
further assess the presence of correlations between me- 
dian activity and age of the PMS regions we have also 
performed two correlation tests (namely the Cox propor- 
tional hazard model test and the generalized Kendall's r 
test), also using ASURV. 

Summarizing the trends for log(Lx) on the left side 
of Fig. |S| as well as the results of the tests, we have 
weak evidence of an increase of log(Lx) in the PMS, for 



|http 



Flaccomio et al.: Activity from the PMS to the ZAMS 



9 



M < O.25M0. More specifically, in the M ^ O.16-O.25M0 
mass bin, p Ophiuchi {\og(Age) = 5.99) may have lower 
levels respect to the outer ONC (log(^5e) = 6.28) and 
the Chamaeleon I region {log{Age) = 6.60). In the follow- 
ing two mass ranges there is no evidence for a temporal 
evolution of median Lx, with the notable exception of 
NGC 2264 (log(^5e) = 6.34), for which the medians, or 
the high luminosity tail of the luminosity functions if the 
median are not defined, are higher than that of the other 
regions of similar age (significantly so in the 0.5-1.0 Mq 
range). Similar lack of evidence for a PMS evolution is 
retrieved for the two highest mass bins, although the two 
population tests indicate a significant lowering of Lx be- 
tween the age of Chamaeleon I and rj Chamaeleontis (the 
two oldest PMS regions) in the 1-2 Mq range. In all mass 
bins the results of the correlation tests do not evidence 
any overall trends of Lx with age in the PMS. A signifi- 
cant decrease of Lx from the PMS to the age of the two 
ZAMS clusters is observed at all masses. 

The trends of log(Lx/Lboi) are somewhat more defi- 
nite. Apart from the lowest mass bin, for which the ev- 
idence from the two population tests is only at the 90% 
level, we observe, for M < IMq, a significant trend of 
increasing log(Lx/Lboi) in the early PMS phase followed 
by a subsequent leveling at the saturation level, ~ —3. 
This is clearly indicated both by the two population tests 
and by the correlation tests (with the exception of the 
0.5 — l.OAf0 for which the result of the correlation tests 
are only at the ^ 85% level, increasing to 97% if we only 
consider the four youngest regions). The saturation level 
is likely retained down to lOOMy for M < 0.5 - I.OMq, 
while at higher masses we observe a decline respect to our 
oldest PMS region. In the two most massive bins all the 
tests do not evidence a statistically significant evolution 
of PMS activity as measured by Lx/Lboi- 

5. Discussion 

We now attempt a physical interpretation of our obser- 
vational results. More specifically we discuss whether the 
observed activity levels can be described with the same 
ideas that are successful for more evolved MS stars, i.e. 
with the mechanism of a stellar dynamo driven by stellar 
rotation plus convection. 

A recent observational study of the dependence of ac- 
tivity on r otation and convection in the MS, has been per- 
formed bv lPizzolato et all l|2003() . who propose two possi- 
ble descriptions of stellar activity levels. 

The first scenario relates the observed X-ray luminos- 
ity (Lx, in ergs/sec) with the stellar rotational period 
(Prot, in days): 

Lx « 10^0 P~2 

This relation is valid, within uncertainties, for rota- 
tional periods such that Lx/Lboi 10^"^ (where Lboi 
depends, on the MS, on stellar mass), while for shorter 
Prot the relation saturates at Lx ~ lO'^^Lboi. Stars of all 



masses therefore follow the same relation of Lx with Prot 
when activity is not saturated but saturate at different 
values of Lx according to their bolometric luminosities. 
In this description convection does not enter explicitly, 
although the presence of convection is necessary for the 
presence of coronal activity, i.e. the relation holds only for 
stars that possess a convective envelope. 

The second description involves the Rossby number, 
i.e. the dimensionless ratio between the stellar rotational 
period and the convective turnover time: Rq = Prot/Tconv 
The observational data, for stars in a fairly wide range of 
masses, can be adequately described by the relation: 

i^L « (2) 

This holds for large enough values of Rq and satu- 
rates, at Lx/L\yo\ 10~^ for smaller Rossby numbers, i.e. 
for short rotational periods or long convective turnover 
times. In this description, like in the previous, stars of all 
mass follow the same relation in the non-saturated regime; 
however, contrary to the previous case, here both the satu- 
ration level and the value of the independent variable, Rq, 
at which saturation occurs are independent from mass. We 
therefore reach an unifying picture at the expense of the 
introduction of t wo additional va r iables , Lboi and Tconv 

As shown by IPizzolato et alJ l)2003|) the two descrip- 
tions are, in the MS phase, equivalent (within observa- 
tional uncertainties): the theoretical stellar structure mod- 
els from which the values of Tconv can be derived indicate 
that, for MS stars in a wide range of mass, the relation 

— 1/2 

Tconv oc L^i^j holds, at least approximately, with a pro- 
portionality constant that is quite independent of stellar 
mass. Equation|21is then approximately equivalent to Eq. 
n By considering only MS stars, it is therefore difficult 
to establish which of the two relations gives a more accu- 
rate description of the data and thus, to shed light on the 
physical mechanisms responsible for activity and for the 
saturation phenomena. 

In the following we compare these two descriptions to 
our PMS data. We will try to decide whether any of the 
two can explain the observations, either entirely or by in- 
troducing so me ad hoc modification . 

Following IPizzolato et al.1 ll2003l). we wil l refer to the 
PMS evolutionary models of I Ventura et al.l l)l998f) : these 
models compute the convective turnover time, whose evo- 
lution in the PMS is essential for our study. Treatment of 
convection is notoriously o ne of the most difficu lt points 
of stellar structure models: IVentura et alJ l)l998|) use two 
alternative convective theories/approximations: Mixing 
Length Theory (MLT) and Full Spectrum of Turbulence 
(FST). For each set of models they compute two differ- 
ently defined values of Tconv : one computed at the bottom 
of the convective zone (not defined for fully convective 
stars), the other throughout the whole convective zone. 
The absolute values of Tconv as a function of stellar mass 
and age depend on its definition and on the choice of con- 
vective theory (MLT or FST). However, the trends with 
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mass and age are similar in all cases. In the following, 
for our argumentation, we will be mainly interested in 
the qualitative trends of Tconv relative to the MS values, 
so that our conclusions do not depend on the particular 
model and/or definition. First of all we note that the in- 
verse relation between Lboi (directly related to mass at 
a given age) and Tconv which holds in the MS does not 
extend unmod i fied to the PMS. Indeed, according to the 
IVent^ira et alJ l|l998() models, moving to younger ages, a 
star of a given mass increases both its Lboi and its Tconv, 
and this latter is more than one order of magnitude larger, 
for IMq stars, at the age of the PMS regions considered 
in this paper, respect to the MS value. The two equations, 
Hand 121 that in the MS are equivalent, will then produce 
in this case different and identifiable predictions. 



5.1. Implications for PMS stars 

We now discuss the implications of the two descriptions of 
activity summarized above for our PMS clusters, with par- 
ticular reference to the ONC case, for which a rich sampl e 
of rotational periods is available fcf. iHerbst et alJ l2002'l. 
We recall that no relation between activ ity and rotation 
has been observed to date iii the ONC ijPeigelson et al.l 
120021 iFlaccomio et all l2003bl) as well as in the other re- 
gions discussed in this work. 

The most direct consequence of the our first scenario 
(cf. Eq. ^ when applied to PMS stars is that, due to the 
larger bolometric luminosity of these latter respect to MS 
stars, the saturation level is predicted to be reached at 
higher values of Lx and smaller rotational periods. Stars 
with measured rotational periods in the ONC sample (Prot 
up to ^ 20 days) should therefore not be saturated. This 
however contrast with the absence of a Lx - Pmt corre- 
lation and with the high, close to saturation, values of 
Lx/Lhoi observed in all our star forming regions and, in 
particular, for stars more massive than O.SMq. For exam- 
ple, excluding the young p Ophiuchi members, more than 
half of the stars in the 0.5 — I.OMq mass range of each of 
our star forming regions have \og{Lx/ L\jo\) > —3.5. Using 
Eq.^and the Lboi values predicted by the SDF evolution- 
ary tracks, we calculate that log(Lx/iboi) > —3.5 would 
imply, for a M — I.OMq star at the age of the ONC, 
rotational periods shorter than 0.8 days. Likewise, satu- 
rated stars (log(Lx/iboi) = —3.0) should have Prot < 0.4. 
However, ^97% and ~99% of ONC stars have measured 
periods longer than 0.8 and 0.4 days respectively. We can 
then exclude that in the PMS Lx is related to Prot by Eq. 

m 

We now turn to Eq. [3 and to the possible role of 
the convective turnov er time, Tconv As already noted the 
IVent^ira et alJ (^2^ models predict that at the age of 
the ONC and younger, the Tconv of solar mass stars is 
more than 1 order of magnitude larger respect to the MS 
value. Assuming that Eq. [21 holds and that saturation oc- 
curs at the same Rq as for MS stars, and recalling that 
MS solar mass stars saturate at log(Prot) ~ 0.2 (with P^t 



measured in days, cf. IPizzolato et aljEoOSj) . an increase 
of an order of magnitude in Tconv translate in an increase 
of the same amount in the rotational periods at which 
saturation occurs. We therefore expect PMS solar mass 
stars with log(Prot) ^ 1-2 to be saturated. Given the sim- 
ilarity of convective turnover times of stars of different 
masses at young ages, saturation periods of the same or- 
der are expected for lower mass stars. Now, log(Prot) — 1-2 
corresponds to periods of about 16 days, i.e. longer than 
typically measured in the ONC. This picture would there- 
fore predict that virtually all PMS stars are saturated. 
Therefore they should not not exhibit a relation of activity 
vs. rotation and they should all have log(Lx/iboi) ~ — 3- 
Such a prediction is almost consistent with the observa- 
tions: we do not indeed observe a correlation of activ- 
ity with rotation and for most star forming regions we 
do observe almost saturated activity levels. This is es- 
pecially true for the older PMS regions (NGC 2264, the 
Chamaeleon I and rj Chamaeleontis associations) and for 
the more massive stars (cf. Fig. Eland IHJ. However, quite 
clearly, the youngest stars, especially at low masses, show 
somewhat lower Px/^boi levels. For the ONC sample of 
rotational periods we have no evidence that rotational pe- 
riods of low mass stars are longer than for more massive 
stars. On the contrary there is likely an opposite depen- 
dence of median P^ot with mass, the median log(Prot) be- 
ing 0.4 (from 15 stars) and ~ 0.9 (from 49 stars) in the 
O.I-O.I6M0 and O.5-I.OM0 mass range, respectively. We 
therefore conclude that, although Eq. (21 i.e. the relation 
between Lx/Lboi and Rq, does better than the relation 
between Lx and Prot (Eq. CJ, it does not by iteeZ/ explain 
the PMS observations. 

One tentative explanation for the low ix/^boi stars 
may invoke the sti ll controyersial existence of the super - 
saturation regime ijRandichL EflOfj jPeigelson et al.[l2002l) . 
There is some evidence of a third regime in the ix/-^boi 
vs. Rossby number relation obtained using MS stars, other 
than the normal inverse relation at large Rq and the satu- 
ration regime: activity of stars with very low Rossby num- 
bers, i.e. very fast rotations or long convective turn-over 
times may show a direct correlation with Rq. The faster 
average rotation of lower mass stars in the ONC, which 
indeed have lower average ix/iboi, seems consistent with 
such an explanation. However, low mass stars in the ONC 
span the whole range of measured rotational periods and 
activity is not correlated, either directlv or inversely, wit h 
rotational period (cf. Fig. 2 in iFlaccomio et al]l2003bl) : 
this holds even when the analysis is restricted to nar- 
row mass ranges in order to reduce any possible difference 
in stellar structure (i.e. Tconv)- Moreover the increase in 
activity in the first few Myr observed at a given mass 
(Sect. 01) would imply, in this hypothesis, a spin-down 
of stars in the early PMS. Although theoretical spec- 
ulations suggest that this could indeed happen due to 
non-instantaneous disk locking (,Hartmann , .2002) . obser- 
vational studies of the rotational ey olution in the ONC 
and NGC 2264 l|R,ebiill et all l2002a^ rather suggest that 
stars keep constant period as they evolve down the convec- 
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tive tracks. Given the lack of firm observational evidence 
of the super-saturation regime, here we do not further ex- 
plore this hypothesis but rather present and discuss an 
alternative explanation. 

5.2. The disk hypothesis 

We now introduce a simple model that tries to account 
for the observed evolution of activity of low mass stars 
in the early PMS phase. Given the persistent lack of ro- 
tational data for the greatest part of our stellar sample, 
the uncertainties in the X-ray data and in the theoretical 
evolutionary models, we keep our model extremely simple 
and, at this point, consider it only as a hint for further 
observational and theoretical work. 

The observational and theoretical bases of our model 
are the foUowings: 

— Most PMS stars, especially in the older regions and/or 
not of very low mass have almost saturated activity 
levels: log{Lx / L^oi) ~ —3.0. A relation between ac- 
tivity and rotation is generally not observed. As dis- 
cussed at the end of the last section, the relation of 
Lx/Lhoi with Rossby number (Eq. [SJ suggests that 
activity should indeed be saturated. 

— Stars with circumstellar accretion disks, in several 
star forming regions, show lower median activity 
levels (both in Lx and L^/ Lhoi) respect to those 
with w eaker or no evidence of a disk and /or ac- 
cretion JStelzer fc Neuhauseil l200lt iFlaccomio et al.L 

The origin of the lower activity level is at 
present unknown. 

We will make the following assumptions: 

1) Departure from saturation is due exclusively to the 

presence of disks, stars that have already shed their disks 

having saturated log(Lx/iboi)- 



scale Tdisk; we can express the fraction of stars that retain 
their disk, /disk, at a time t by: 



log(Lx/ibol)sat. 



(3) 



2) Stars with a disk all have the same value of 
log(Lx/iboi) which we will take to be -5.0. 



H, 



disk 



log(Lx/ibol)disk 



(4) 



Although this value is quite arbitrary, it is compat- 
ible with the lowest values of log(Lx/iboi) of low mass 
(< 3Mq) ONC and p Ophiuchi stars. According to max- 
imum likelihood distributions, only ^2 0% of both the 
high a ccretion ONC stars (as defined bv IFlaccomio et aD 
l2003bt) and of all p Ophiuchi members in our sample, have 
log(Lx/iboi) below this threshold. Note also that in a 
more realistic model which allows for a continuous range 
of ix/iboi values and of disk evolutionary stages, i/disk 
would correspond to log(Lx/iboi) of stars with the least 
evolved accretion disks. 

3) The dissipation of circumstellar disks can be de- 
scribed by a stochastic process with characteristic time- 



J disk — c 



(5) 



We can then compute the mean log(Lx/iboi) as a 
function oi t, < H{t) >. 



< H{t) i/sat. + (i^disk - ifsat.) e-*/"^-^ 



(6) 



This model has two unknown parameters: iJdisk, which 
we have fixed at -5, and Tdisk which we keep as a free pa- 
rameter. Hsat is here considered as fixed by the observa- 
tions at -3, roughly corresponding to the highest values 
observed in the PMS and equal to the saturation level ob- 
tained for fast rotating MS stars. The right column of 
Fig. IHl shows, superposed on the already described ob- 
served trends of log(-Lx/iboi) with age for different mass 
ranges, three realization of the model corresponding to 
three values of the disk dissipation time-scale, 1.0, 1.5 and 
2.0 Myrs. We note two facts: 

1) The agreement between data for the PMS regions 
and models is qualitatively good, also considering the ex- 
treme simplicity of our assumptions. The models do not, 
and are not supposed to, reproduce the fall in Lx/^boi 
on the MS, which can be attributed either to wind-driven 
spin down of the stars or, for the 2 — ZMq bin, to the loss 
of the convective region. 

2) The observed evolution of Lx/-^boi is somewhat 
faster for more massive stars. The model seems to suggest 
that Tdisk depends on stellar mass, ranging from 2.0 
Myrs for M - 0.21/© to - 1.0 Myrs for M - 0.75Mq, 
where the quoted masses are simply the centers of the 
second and fourth mass bins. 

We therefore seem to have derived an indication on the 
circumstellar disk lifetime. How do these results compare 
with totally independent results obtained from specific 
studies performed at infrared wavelengths? iHaisch et alJ 
analyze a sample of six clusters and determine, 
from a linear fit of disk fraction vs. cluster age, a disk 
lifetime of about 6Myr, this being the age at which essen- 
tially all the stars loose their near-lR detectable disks. A 
similar analysis performed by Hillenbrand (,2Q03i) using a 
larger sample of star forming regions, suggests instead an 
inner disk lifetime of only 2-3 Myr. An independent indi- 
cation on the timescale f or the evolution of in ner disks has 
been recently derived bv lHerbst et alJ (2Qq3) from the ro- 
tational periods of ONC members. Interpreting the period 
distribution in the context of the disk locking hypothesis, 
they conclude that the timescale for ONC stars to be re- 
leased from their disk locking is 0.5 — 1.0 times the ONC 
age. If the two median stellar ages reported in Fig. ^ for 
our ONC sub- regions are used this translates to 0.7 — 2.0 
Myr. 

Our best guess for Tdisk derived from X-ray data, i.e. 
l-2Myrs., falls within the range of the above estimates. 
We note that all these estimates are on one hand sub- 
ject to significant uncertainties, and on the other hand 
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may trace different, though related phenomena such as the 
presence of the inner disk, disk locking and mass accretion. 
The order of magnitude agreement is therefore suggestive. 
Although we have identified the evolutionary time-scale of 
Lx/^boi with the disk dissipation time-scale, our model is 
actually rather general: it only assumes a two level popu- 
lation and a fixed transition probability per unit time from 
the low to the high level. We have tentatively identified the 
low level with the pres ence of disks because recent studies 



l|Flaccomio et al.L l2003bna) indicate that stars surrounded 
by accretion disks have lower levels of log(Lx/-^boi)- If 
the lowered activity level is in fact due to accretion it- 
self as opposed to the presence of a disk, the time-scale 
for the evolution of activity might correspond to the time- 
scale for the evolution of accretion, which might be shorter 
than the time-scale for the disappearance of IR detectable 
disks. The rather good agreement of our X-ray derived 
time scale w i th the disk locking times scale reported by 
iHerhst et ai .1 (|2nn2r) is particularly suggestive as disk lock- 
ing implies a modification of the magnetic field geometry, 
conceivably affecting the X-ray emitting plasma. 

As for our second result, i.e. the dependence of the 
X-ray derived Tdisk from stellar mass, we note that this 
dependence is also suggested, and i n the same direction 
found here, bv several IR studies fe.g. lHaisch et alJboOlb! : 
iHillenbrand et al]ll998l) . thus giving further support to 
our interpretation. 

More X-ray and optical/IR data are needed to fully ex- 
plore the relation between magnetic activity and accretion 
and/or disks. 

6. Summary and Conclusions 

We have studied the evolution of X-ray activity in the first 
stages of stellar evolution, from the early PMS (~ 1 Myr.) 
to the ZAMS (~ 100 Myr), using published and archival 
data on five star forming regions and two young ZAMS 
clusters. In all cases we have reviewed the X-ray and op- 
tical data used to characterize members of the respective 
regions. If necessary, we have repeated the data analysis 
in order to homogenize and update, to the best of our 
knowledge, the assumptions used to derive activity indi- 
cators {Lx and Lx/-^boi) and stellar parameters (masses 
and ages). 

In the early PMS phase we observe that low mass stars 
that posses a convective region show high levels of activ- 
ity, often close to the saturation level (log Lx/-^boi ^ —3) 
found for fast rotating MS stars. We compare these activ- 
ity levels with those predicted by two physically distinct 
scenarios: 

#1) the total energy output of the corona (Lx) is de- 
termined exclusively by stellar rotation (Prot) and satu- 
rates when Lx reaches 10"'^ times the total stellar lumi- 
nosity. 

#2) the fraction of energy output by the corona re- 
spect to the whole energy produced by the star is deter- 
mined by the dimensionless ratio between the characteris- 
tic time of rotation and the characteristic time of convec- 



tion (i.e. the Rossby number, Rq). This relation saturates 
at a fixed value of Rq, corresponding to Lx/Lboi ^ 10^^. 

While, for MS stars, these two scenarios are hardly dis- 
tinguishable, both giving a fair description of the observa- 
tions, they predict distinctively different results if applied 
to 1 — 5 Myr old stars. Scenario #1, in order to explain 
the high levels of activity observed, would require very fast 
rotations, which are simply not observed. Scenario #2 on 
the other hand predicts that stars in the early PMS with 
L'rot ^ 16 d, i.e. consistent with the rotations measured in 
the ONC, should all be saturated. 

Our data seem therefore to favor scenario #2, although 
it does not fully account for the observations: median 
Lx/Lboi appears to be somewhat lower than the satu- 
ration value for young ages (< 5Myr) and low masses 
(< IMq), but tends to reach that level at older ages and 
larger masses. Saturation is then retained up to the age of 
the Pleiades and NGC 2516 (100- 140Myrs) for stars with 
M < 0.5Mq, while activity in more massive stars seems 
to desaturate earlier. We can tentatively explain this lat- 
ter decay of activity from the PMS to the MS with the 
decrease of the convective turnover time (or, equivalently, 
convective zone depth), which is indeed more pronounced 
for more massive stars l|Ventura et al.L Il99i^ and, possi- 
bly, with the angular momentum loss induced by magnetic 
winds (cf. Bouvicr ct a l.,1997.) . The increase of Lx/Lboi in 
the early PMS however is not so readily explained by the 
foreseen evolution of the Rossby number. Supported by 
observational evidence, we put forth a simple model that 
assumes the existence of two populations: one in which 
activity is at the saturation level, the other with lower 
observed activity, which we tentatively identify with stars 
surrounded by circumstellar disks and/or undergoing ac- 
cretion. We then postulate that the low activity popula- 
tion makes a transition to the high activity population 
(e.g. dissipates disks or stops accreting) with a character- 
istics time, T. Observed median activity levels are roughly 
consistent with such a model if r ~ 1 — 2 Myr, being at the 
short end of the interval for higher mass stars (~ IMq) 
and at the long end for the lowest mass ones (~ O.2M0). 
These characteristic times are similar to disk dissipation 
times derived by IR studies, to the timescale for disk lock- 
ing, as derived by modeling the rotational periods distri- 
bution of ONC stars, and to the expected timescale for 
the evolution of circumstellar accretion. Given the consid- 
erable uncertainties in all of these estimates, we restrain 
from any more precise identification of the physical pro- 
cess that determines the evolution of activity. We also note 
that the observed trend of r with mass is consistent, ac- 
cording to recent IR studies, with the disk (and possibly 
with the accretion and disk locking) hypothesis. 
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Fig. 8. Time evohition of median Lx and ix/^boi (left and right respectively). Medians are represented by small 
horizontal segments, with hatched boxes representing uncertainties. Down-pointing arrows indicate upper limits on 
the median and on the lower error (thick and solid arrows respectively). In some cases medians are substituted with 
averages (either defined or upper limits), indicated by error bars. Numeric values above abscissa indicate numbers of 
detections ("Det.") and upper limits ("Upp.") in each of the bin. The three curve in the right hand panels refer to 
three realizations of the simple model described in Sect. |S1 characterized by the disk dissipation time, r<jisk, whose 
value is reported beside the curves in units of 10® years. 



